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Abstract 

The stability of the bound states of the magnetic flux in a Josephson 
resistive lattices is investigated numerically. It is shown that for a simple 
relationship between the geometrical parameters of the lattice the range of 
bias current is the widest. 
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I. Formulation of the problem 



In the present paper we study numerically stationary spatially periodic 
states of the magnetic flux in one- dimensional Josephson junctions, whose 
dielectric layer represents a lattice of resistive inhomogeneities of finite size.0'i 

The mathematical model of such junctions is described by the perturbed 
stationary sin-Gordon equation 

- + JD(x)sin0 + 7 = 0, xe{-R,R), (1) 

supplemented with the periodic boundary conditions 

0(-i?)=0(i?), 0,(-i?)=4(i?). (2) 

Here joix) is given continuous periodic function which models the distribu- 
tion of the amplitude of the Josephson current along the junction; 2R is a 
given length of a spatial wave; 7 is the bias current. 

Let us consider in more details the choice of the function jnix). For 
the homogeneous junctions we have joix) = 1. In an inhomogeneous case 
the function jnix) is defined by the number of the inhomogeneities and 
their shape. The (5-function model of inhomogeneities is the most frequently 
used. In this paper we use inhomogeneities of the shape of equilat- 

eral trapezium (see fig. 1). Such trapezium is characterized by the width of 
its bottom /i and its top a = mfi. In all cases the height of the trapezium 
equals to a dimensionless unitJl 

For m = 1 the inhomogeneity is a rectangular with base /i. The more 

realistic model of an inhomogeneity is one with m > 1. The decay of the 

amplitude of the Josephson current from to 1 in such an inhomogeneity 

Tn — 1 

takes place in the interval of length 6 = — - — n for the particular choice of 
the width /i of the inhomogeneity. 

Most of the calculations in this work were carried out for m = 1.5, i.e. 
for 6 = 1/4. We consider an "infinite" lattice consisting of the described 
inhomogeneities with separation A between their center positions. The sim- 
ple relationship /i < ^ between the values of A, /i and m ensures that two 
adjacent inhomogeneities do not overlap. 

The lattice is symmetric when A — o" = /x; this yields /i = . 

1 + m 
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To examine the stability of a certain solution of nonlinear boundary 
value problem (BVP) (|l|), (0) we solve the Sturm-Liouville problem (SLP) 



V',(±i?) = 0, (4) 
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if) [x] dx = 1, (5) 



where the potential 

q{x) = joix) cos (j){x) (6) 

is generated by this particular solution. The solution is stable if the 
minimal eigenvalue Aminof the SLP with the corresponding potential q{x) 
given by (|^), is positive.u 

To solve numerically the nonlinear problems (0), |^ and |^ - we use con- 
tinuous analog of the Newton's methodS with an optimal step.0 At each 
iteration we discretize the corresponding linearized BVP using the spline- 
collocation second order scheme over a nonuniform grid with condensed spac- 
ing in the vicinity of the inhomogeneities. The detailed description of the 
algorithm and the numerical scheme is given in.§ 



II. Discussion of the numerical results 

Each solution (j){x, A, ^, Nj, 7) of BVP (|l|), (|) depends on the geometri- 
cal parameters A, n and Nj, and on the physical parameter 7 as well. Here 
Nj is the number of inhomogeneities, enveloped by a single wave of length 
2R = NjA. Such waves exist in a real finite size junctions having Nj inho- 
mogeneities symmetric with respect to its center. Here 2R is the length of 
the junction. We note that in this case the eq. (0) is supplemented by the 
boundary conditions of Neuman type 

0,(±i?) = hn. (7) 

Here /i^ is the magnetic field at the boundaries of the junction. 

The potential q{x) of the SLP @ - (|^) is explicitly related to the periodic 
solution of the BVP ([1|), (0) by eq. (|^), and, consequently, the corresponding 
eigenvalues and eigenfunctions depend on the parameters A, /i and Nj. Thus 
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the stability of the bound states of the magnetic flux in the junction also 
depends on these parameters. 

Note that when the bias current 7 = the eqs. (|l|), @ have "vacuum" 
(Meissner's) solutions of the form 

(j){x)=kTT, A; = 0,±1,±2,..., 

for every Nj. The solutions corresponding to even k are stable while the 
solutions corresponding to odd k are unstable. 

Let us consider the influence of the geometric parameters on the possible 
bound states of the magnetic flux in the lattice. 

For a junction with Nj = 1, 2R = A < 27r and 7 = we find just 
two periodic solutions — stable and unstable Meissner's solutions. Such 
junction is "short" for periodic solutions, however, we should note that there 
are nontrivial stable aperiodic solutions corresponding to the eq. (|I]). The 
derivative (x) of one of these solutions is marked by "o" in fig. 2. For 
large enough values of the distance A, nontrivial unstable periodic solutions 
exist in the junction. The most interesting among them is the one, whose 
derivative is marked by "A" in fig. 2. We will conditionally call this solution 
soliton. It remains unstable in a very wide range of the values A and /i; this 
is demonstrated in fig. 3 and fig. 4. The second unstable solution, marked by 
"□" in fig. 2, exists in a real junction with one inhomogeneity only for very 
large values of Hb at the boundaries. 

There is a periodic stable solution in the lattice in the case Nj = 2. Fur- 
ther we refer to this solution as "main periodic soliton". This solution is 
marked by "o" in fig. 2 and fig. 5. If the two symmetric inhomogeneities of a 
junction are separated by large enough distance 6, then the periodic magnetic 
flux (^{x) is a result of the nonlinear interaction between two aperiodic sta- 
ble pinned at the inhomogeneities structures - the fluxon and the antifluxon. 
Similarly, for the case Nj = 1 the unstable solutions are the result of "in- 
teractions" between two unstable fluxons, alienated from the inhomogeneity. 
For these reasons one can make the conjecture that for the case Nj = 1 there 
aren't any nontrivial periodic stable solutions.^ 

The solution, marked by "V" in fig. 5 (the "secondary" periodic soliton) 
is an example of weakly stable bound state, whose energy is approximately 
four times greater then the energy of the main periodic soliton. The stable 
soliton with small Amin exists only for very large values of the parameters 
A and fi of the lattice. Both solitons have symmetric reflections about the 
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Further we refer to them as main periodic antisohton and secondary 
periodic antisohton. The remaining two solutions in fig. 5, marked by "□" 
and "o" are unstable. 

Fig. 6 shows the minimal eigenvalue Amin, corresponding to two stable 
periodic solutions, as a function of the parameter fi when Nj = 2 and bias 
current 7 = 0. Different curves correspond to different values of the distance 
A between the inhomogeneities. It is very important to note, that all the 
curves have maximum in a certain point /i^, dependent on the parameter A. 
This leads that in a lattice with inhomogeneities separated by the distance 
A with width /i^, the time-dependent perturbations of the magnetic flux 
will decay smoothly .0 The location of the points of maximum /i^ for all the 
solutions depends on the value of A and for the main soliton (as for main 
antisohton) is equal approximately to 1.3. 

Apart from the periodic stable solutions in a real junction with two inho- 
mogeneities there are stable aperiodic bound states of the magnetic flux when 
the boundary magnetic fleld hs = 0; this case corresponds to the boundary 
conditions (|^). 

The distribution of the aperiodic magnetic fleld 4>^{x) along the junction 
is shown in flg. 7 (the symmetric reflections are not shown). Here the values 
of the parameters are R = 6, Ni = 2, Hb = 0, 'j = 0. The marked by 
"o" solution is periodic, others are not. The energy of the "complex" bound 
state (periodic soliton), pinned at two inhomogeneities, is much greater than 
the energy of the aperiodic solitons. The minimal eigenvalue of the periodic 
soliton is smaller than the minimal eigenvalues corresponding to the aperiodic 
bound states. 

The minimal eigenvalue Amin as a function of the parameter /i has an 
maximum for the aperiodic bound states as well as for the aperiodic ones. 
This is demonstrated in flg. 8. 

For comparison, the same curve (marked by "o"), corresponding to the 
stable fluxon in the junction having one inhomogeneity, is shown in flg. 9. 

The curves Amin(/^), corresponding to the unstable periodic solutions, are 
shown in flg. 9. Here the parameters of the lattice are Nj, n = 2.4 and 
the bias current 7 = 0. The stable main periodic soliton transforms into 
a unstable solution at the point of bifurcation Bq (the unstable soliton is 
marked by "□" in flg. 5). The secondary periodic soliton is unstable in the 
lattice with so "narrow" inhomogeneities. The secondary soliton transforms 
into another unstable periodic soliton at the point of bifurcation Bi when 
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the second eigenvalue of the eqs. - is also equal to zero. The latter 
unstable soliton is marked by "o" in fig. 5. Increasing of the width of the 
inhomogeneities "lifts" the curve Amin(/^) upwards and, if fi is greater than a 
certain critical value, the secondary periodic soliton becomes stable. 

Let us consider the dependence of the bound states of the magnetic fiux 
on the bias current 7. It is well-known that in the inhomogeneous junction 
the stable solutions lose their stability when the absolute value of the bias 
current is increased. In the periodic case the process of "destruction" of 
the magnetic field by the bias current is demonstrated in fig. 10. It can be 
seenjli that the positive values of the current 7 "compress" the field in the 
center between two inhomogeneities, while the negative values of 7 push it 
to the adjacent inhomogeneities or, in the case of a real junction with two 
inhomogeneities, to the boundaries of the junction. 

The curves Amin(7) for all found stable solutions as well as unstable, in 
the case Nj = 2, A = 6, fi = 2.4 are demonstrated in fig. 11. The curves, 
corresponding to the stable solutions, have their symmetric refiections about 
the fine 7 = 0. As before, Bq denotes the bifurcation points, corresponding 
to the roots of the equation Amin(7) = 0. 

Let 7„ and 7j, are the roots of this equation (the points Soa and -B06 
correspond to this roots in fig. 11 for the main soliton/antisoliton). Let us 
designate the difference 

A7 = 7f, - la 

interval of stability of the main periodic soliton with respect to the bias 
current 7. Generally, the value of A7 is a function of the parameters A 
and /i. The interval of stability A7 for the main soliton as a function of 
the parameter /i, when the value of the distance A is fixed, is shown in the 
fig. 12. These curves can be approximated by the cubic polynomials with 
coefficients, depending only on A, in the whole domain where /i is defined 
physically. For a given A the points of maximum determine the width of 
the inhomogeneities, so, that the soliton has the maximal interval of stability 
A7 with respect to the current 7. 

Note, that the shape of the curve /i^(A) can be fitted by a fine /is(A) = 
a A + b with a good accuracy. This is the line, marked by "o" in fig. 13. The 
approximate values of the parameters of this line, calculated numerically, are: 
a ^ 0.254, b ^ -8.5. lO^^. 

Having in mind the above results we can construct an "optimal" junction 
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with geometric parameters related by 

In this case we have Nj — 2, Aopt ~ 4//c- Thus, the "optimal" distance 
between the inhomogeneities is Aopt ~ 5.2. 

III. Conclusions 

We show numerically that in the Josephson junctions with resistive inho- 
mogeneities the stable periodic bound states of the magnetic flux, pinned on 
the inhomogeneities, have a certain minimal wavelength. The dependence 
of the stable bound states on the width of the inhomogeneities and the 
distance A between them is studied. In particular, it is demonstrated that 
the curve XminilJ') for a fixed A has maximum /j,^, i.e. in such a junction the 
time-dependent perturbations of the magnetic flux decays smoothly. The 
dependence of the interval of stability A7 on the parameter /j, is also studied. 
For a certain relationship between the with of the inhomogeneities and their 
separation ("optimal" lattices) the interval of stability A7 is the widest. 
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V. Figure captions 



1. Geometrical model of a resistive inhomogeneity. 

2. Some distributions of the magnetic field 0a. (x) along the junction. The 
solution, marked by "o", is stable aperiodic soliton of the magnetic field in 
the case Nj = 1; the solutions, marked by "□" and "A" are unstable periodic 
(case Nj = 1); marked by "o" solution is stable periodic (case Nj = 2). 

3. Minimal Eigenvalues Xmin of the SLP (^) - (^ versus the distance A 
between the inhomogeneities for the solutions, marked by "□" and "A" in 
fig. 2. 

4. Minimal Eigenvalues Amin of the SLP (^ - (^ versus the width fj, of 
the inhomogeneities for the solutions, marked by "□" and "A" in fig. 2. 

5. Stable (marked by "o" and "□") and unstable (marked by "o" and 
"□") distributions of the magnetic field 4>x{x) along the junction with pa- 
rameters Nj = 2, A = 9.5, /i = 2.4. Here 7 = 0. 

6. Minimal eigenvalue Amin of the SLP (^) - (||) versus the width /i of the 
inhomogeneities for the stable solutions from fig. 5. Here the parameter A is 
fixed and the point /i^ corresponds to the "optimal" width. 

7. Stable bound states of the magnetic field 0^(x) corresponding to the 
aperiodic boundary conditions (|^). Here Nj = 2, A = 6, fi = 1, •y = 0, 
hs = 0. The quantity Acj) = (j){R) — (f){—R) is the full magnetic flux through 
the junction. 

8. Minimal eigenvalue Amin of the SLP (H) - (|5D versus the width /i of the 
inhomogeneities for the stable solutions from fig. 5. Here the parameter A is 
fixed. 

9. Minimal eigenvalue Amin of the SLP (|^) - (^ versus the distance A 
between the inhomogeneities. Bq and Bi are the points of bifurcation. 

10. Variation of the shape of the stable soliton as the bias current is 
increased. Here A?^/ = 2, A = 6, = 2.4. 

11. Minimal eigenvalue Amin of the SLP @) - (||) versus the bias current 
7. -Boa and -Bob — for the periodic soliton from fig. 7. 

12. Interval of stability A7 = 7;, — 7^ versus the width /i of the inhomo- 
geneities. The parameter A is fixed. 

13. The width /i^,, corresponding to the largest interval of stability versus 
the distance A in the cases Nj = 2 and Nj = 3. 



8 



References 

^N.V. Alexeeva, T.L. Boyadjiev, Bulgarian Journal of Physics, v. 24, 1997, 
No 1-2, p. 63. 

2T.L. Boyadjiev, N.V. Alexeeva, Comm. JINR, P17-97-28, Dubna, 1997 (In 
Russian) . 

^Yu. S. Gal'pern, A.T. Filippov, JETF (Soviet Journ. Teor. Exper. Fiz), v. 
86, No 4, 1984, p. 1527 (In Russian). 

''I.L. Serpuchenko, A.V. Ustinov, Solid State Comm., 68, 7, 1988, p. 693. 

^B.A. Malomcd, A.V. Ustinov, Fizika Nizkih Temp., 15, (1989), p. 1128. 

^B.A. Malomed, A.V. Ustinov, Phys. Rev. B, v. 49, No 18 (1994), p. 13024. 

^R. Fehrenbacher, V.B. Geshkenbein, G. Blatter, Phys. Rev. B, v. 45, No 
10 (1992), p. 5450. 

«B.H. Larsen, J. Mygind, A.V. Ustinov, Phys. Lett. A, 193 (1994), p. 359. 

^A. V. Ustinov, Josephson vortices in distributed superconducting structures, 
D.Sc. Thesis, Tchernogolovka, 1994 (In Russian). 

^°E.P. Jidkov, G.I. Makarenko, I.V. Puzynin, In "Physics of elementary par- 
ticles and atomic nucleus", Moscow, Atomizdat, 1973, v. 4, No 1, p. 127 
(In Russian). 

^^V.V. Ermakov, N. N. Kalitkin, JVMiMF, v. 21, No 2, 1981 (In Russian). 

^^A.V. Filippov, 1992 (private communication) . 

^^D.W. McLaughlin, A.C. Scott, Phys. Rev. A, 18 (1978), p. 1652. 



9 



